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When we look at a non-deterministic system, it is important to check the variance and covari- 
ance of the output y. It determines how certain we are of our results. When we look at the 
output, there are two types of question we might encounter. The type is which the transient 
state (when the system just starts ) matters. The second type we only look at the steady state. 
This assumes that we can let the system warm up for some time t. In this paper we are trying 
to find the mean/covariance of a system at the steady state. 

Let us consider this system: 



where ui and v are unit strength white noises, x and v are independent. C = l/10w„ = 2. 

Let's first find the mean of x and y 

we know that x can be solved as 




y = ( 1 1 )x + v 



x(t)= f_ x e A ^- T ^B^(T)dT VteJR 



the expectation of x(t) would then be 



E[x(t)]=E^J t _ cc e A{t - T ">BLo( T )dT Wen] 

e A(t ~ T) B E[lo(t)] dr \/teH 




now if we want to find the expectation of y 



y = ( i i )x + v 



E[y] = E[( l l )x + v] 
E[y] = (i i)E[x]+E[v] 
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this is also pretty simple since we already calculated that the expectation of x = 0 and v 
is white noise so E[v] = 0 



E[y} = 0 

Next we want to find the covariance of x and y 

Cov[ x ] = E[ z(f)x(i)* ] = P 
P is the solution to the lyapunov equation such that 
AP + PA* + BWB* = 0 

You can solve for P by hand, but I do not recommend it instead try these code in 

matlab: 

z = 1/10.0: 

w = 2; 

A = [0 l;-w2 -2*z*w]; 

Q = [1 i;i l]; 

lya = lyap(A,Q) 

And we get the solution P as 

lya = 

1.8625 -0.5000 
-0.5000 6.2500 



Next we want to find the Covariance of y 

Cov[i/] = E[ -E[i/(t)]E[ y(t)}* 

since the expectation of y is zero we only have to worry about the first term. 
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Cav[y] = E[ i/(t)i/(t)*] 

E[ i/(t)i/*(t)] = £?[(( 1 1 )! + «)(!*( J ) +v*) ] 

E[ y(t)y*(t)}= (i i )E[a:a:*]( J ) + J E[w*a:*]( \ ) + ( i i )B[a;t;*] + ] 
Previously we stated that v and x are independent so their covariance would be zero 
E[y(t)y*(t)}= (i i)E[xx*](\)+l 

e[ „(*)»•(*)]=(! --)(;) + i 

if you add these up in matlab you would get 

Cy = 

8.1125 



Now let's try to find the expectation on this 
E [■ lim I £ x(t)'Qx(t)dt] 

IT— >oo J 

where Q 

«=(?!) 

using this property 

E [- lim \ J 0 T z(t)'z(t)dt^ = trace(P z ) 

z in this case would be Q 1/2 x(t) 

z=Q^x( t )=(i iy/Mt) 

P^ElQWxWx'WQ* 1 '*] 
Y z = Q 1 l 2 E[x{t)x*{t)]Q* 1 l 2 

if we put this into matlab we would get 
trace(P z ) = 15.225 
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